Draft version December 12, 2008 

Preprint typeset using I£T^X style cmulateapj v. 10/09/06 



FASHIONABLY LATE? BUILDING UP THE MILKY WAY'S INNER HALO 

Heather L. Morrison 

Department of Astronomy, Case Western Reserve University, Cleveland, OH 44106 

Amina Helmi 

Kapteyn Astronomical Institute, University of Groningen, PO Box 800, 9700 AV Groningen, the Netherlands 

Jiayang Sun, Peng Liu and Rongfang Gu 

Department of Statistics, Case Western Reserve University, Cleveland, OH, 44106 

John E. Norris 

Research School of Astronomy and Astrophysics, The Australian National University, Mt Stromlo Observatory, Cotter Road, Weston 

ACT 2611, Australia 

Paul Harding 

Department of Astronomy, Case Western Reserve University, Cleveland, OH 44106 

T.D. KlNMAN 

National Optical Astronomy Observatories, PO Box 26732, Tucson AZ 85726 

Amanda A. Kepley* 

Department of Astronomy, University of Wisconsin-Madison, 475 North Charter Street, Madison, WI 53706 

Kenneth C. Freeman 

Research School of Astronomy and Astrophysics, The Australian National University, Mt Stromlo Observatory, Cotter Road, Weston 

ACT 2611, Australia 

Mary Williams 

Research School of Astronomy and Astrophysics, The Australian National University, Mt Stromlo Observatory, Cotter Road, Weston 
ACT 2611, Australia and Astrophysikalisches Institut Postdam, An der Sternwarte 16, D-14482 Potsdam, Germany 

Jeffrey Van Duyne 

Department of Astronomy, Yale University, P.O. Box 208121, New Haven, CT 06520 
Draft version December 12, 2008 

ABSTRACT 

Using a sample of 246 metal-poor stars (RR Lyraes, red giants and RHB stars) which is remarkable 
for the accuracy of its 6-D kinematical data, we find, by examining the distribution of stellar orbital 
angular momenta, a new component for the local halo which has an axial ratio c/a ~ 0.2, a similar 
flattening to the thick disk. It has a small prograde rotation but is supported by velocity anisotropy, 
and contains more intermediate-metallicity stars (with —1.5 < [Fe/H] < —1.0) than the rest of our 
sample. We suggest that this component was formed quite late, during or after the formation of the 
disk. It formed either from the gas that was accreted by the last major mergers experienced by the 
Galaxy, or by dynamical friction of massive infalling satellite(s) with the halo and possibly the stellar 
disk or thick disk. 

The remainder of the halo stars in our sample, which arc less closely confined to the disk plane, 
exhibit a clumpy distribution in energy and angular momentum, suggesting that the early, chaotic 
conditions under which the inner halo formed were not violent enough to erase the record of their 
origins. The clumpy structure suggests that a relatively small number of progenitors were responsible 
for building up the inner halo, in line with theoretical expectations. 

We find a difference in mean binding energy between the RR Lyrae variables and the red giants 
in our sample, suggesting that more of the RR Lyraes in the sampl e belong to the outer halo, and 
that the outer halo may be somewhat younger, as first suggested bv lSearle fc Zinnl (|1978l ). We also 
find that the RR L yrae mean rotation i s more negative than the red giants, which is consistent with 
the recent result of ICarollo et al.l (|2007f ) that the outer halo has a retrograde rotation a nd with the 
difference in kinematics seen between RR Lyraes and BHB stars bv lKinman et al.l (|2007T) . 
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Subject headings: Galaxy: halo — Galaxy: formation — Galaxy: evolution — Galaxy: kinematics 
and dynamics 



1. INTRODUCTION: RECONSTRUCTION OF HALO 
HISTORY 

One of the central aims of the study of old popula- 
tions in the Galaxy is the reconstruction of its history 
from their properties today. We will begin here by con- 
sidering the differ ent approaches to mo deling the struc- 
ture of the halo. iSearle fc Zinnl (|1978t ) took an impor- 
tant first step toward our current understanding of the 
hierarchial build-up of galaxies when they noted that 
"the loosely bound clusters of the outer halo have a 
broader range of age than the more tightly bound clus- 
ters." Perhaps because of our taxonomic heritage, sub- 
sequent stellar populations workers have tended to de- 
scribe the complex structure of the stellar halo in terms 
of two components. For example, in one of the first at- 
tempts at "Armchair Cartography" (the use of a local 
sample of halo stars to reconstruct the h alo density dis- 
tribution) ISommer^ Larsen fc Zhenl (|1990f ) inferred a two- 
component halo, with one component quite flattened. 
However, the study was limited by the small sample and 
the relatively inaccurate distances and velocities then 
available. 

More global studies of the halo using kinematically 
unbiased samples also found v a riants on this two- 
compo nent halo. Kin man et alj (|1965f ) and lHartwickl 
(1987) used RR Lyrae variables to show that the flat- 
tened component dominates in the inner halo, and the 
spheri cal c omponent in the out er halo. iPreston et al.l 
(|1991[ ) and iKinman et all (|1994f ) found a similar spatial 
distribu tion using blue horizo ntal branch stars and RR 
Lyraes. IChiba fc Beers! (|2000t ) used a larger, local sam- 
ple and confirmed the result s of previous studies. Most 
recently. ICarollo et all (|2007h used more than 20,000 lo- 
cal stars from SDSS to further quantify this "halo di- 
chotomy" . They find that the outer halo is significantly 
more metal-poor than the inner halo, and has a retro- 
grade mean rotation. 

The influenti al work of IZinnl (119931) (building on the 
earlier work of lSearle fc Zinnl Il978h developed the two- 
component halo idea further, using globular cluster kine- 
matics. Zinn subdivided the halo globular clusters into 
different groups using horizontal branch morphology (as 
a proxy for age) and found different spatial distributions 
and kinematics for the two groups. The Old Halo clusters 
in the inner halo show a flattened distribution and have 
prograde rotation (although much less than the disk's). 
The Younger Halo is spherical and has no net rotation. 

With the excep tion of th e work of 

ISommer-Larsen fc Zhenl (|1990l ) and IChiba fc Beersl 
1)20001 ) . all of the works which derived a flattened inner 
halo were based o n samples of s t ars m ore distant than 
~1 kpc. However. IPreston et alj (|1991l ) noted a curious 
anomaly: when they used their density model (including 
an inner halo with c/a=0.7 at the Sun) to predict the 
density of RR Lyraes in the solar neighborhood, the 
model predictions fell seriously short: by a factor of at 

* Current address: Department of Astronomy, University of Vir- 
ginia, Charlottesville, VA 22904-4325 



least two, even though the local sample is unlikely to be 
complete at very low latitudes. Preston et al. followed 
this fact to its logical conclusion and suggested that the 
halo has an additional highly flattened component found 
only very close to the plane. This paper will confirm the 
detection of the highly flattened halo component with a 
new sample of halo stars which have a median distance 
of 1 kpc. 

How can we explain these results in the context of cur- 
rent models of galaxy formation? There have been sig- 
nificant advances in recent years which have developed a 
more nuanced picture of galaxy formation in a number 
of areas. The concept which has changed least over the 
years is that of disk formation: a smooth, dissipational 
collapse which conserves the angular momentum of the 
infalling material seems the only way to produce large 
galaxy disks like our own. Halo formation was likely 
much messier. Our concepts here have changed more 
with time, as the hierarchical paradigm for structure for- 
mation has permeated the field: we now understand that 
accretion is fundamental to the formation of the dark 
halo. Classical expectations are that this accretion is 
violent in its early stages: the rapid variations in the 
gravitational potential caused by early mergers should 
smooth out any record of the discrete origins of early 
halo stars. We will see below that this is not entirely 
correct. 

Before the disk forms, halo stars could be produced by 
th e mergers of gas-rich o bjects (the "dissipative mergers" 
of iBekki fc ChiblJl200lh . (Only minor accretions can oc- 
cur after the present-day disk forms, because otherwise 
the fragile kinematically cold disk would be destroyed). 
The stellar halo can also be built up by accretions of 
larger objects, helped by dynamical friction (which was 
likely the way the progenitor of the globular cluster us 
Centauri arrived in the inner halo) , or by the simpler ac- 
cretion of smaller satellites which did not feel dynamical 
friction. We will see that there are opportunities for all 
of these processes to build up a flattened inner halo. 

The use of integrals such as energy and angu- 
lar momentum, which are conserved under some con- 
ditions as the Galaxy evolves, can be a powerful 
tool for tracing the Galaxy's his tory. For example, 
iLvnden-Bell fc Lvnden-Beill (| 19951) and Uohnston et ail 
(|1996| ) showed how energy can be used to i dentif y 
stars with commo n orig in s, and iHelmi et al.1 lll999T ) , 
I Helmi fc de Zeeuwl (pOOOj) . IChiba fc Beersl (|2000l ) and 
IHelmi et al.l (|2006) considered d ifferent components of 
angular momentum and energy. IHelmi et all (|1999h il- 
lustrate the different information available from study 
of stellar velocities and of angular momenta: their Fig- 
ure 2 shows that the debris from a disrupting satellite 
is all found in the same place in the angular momentum 
diagram, but that different streams from this satellite 
may be seen with distinct velocities (moving towards and 
away from the plane, for example). 

In this paper we analyse a sample of metal-poor stars 
from the solar neighborhood which is remarkable for the 
accuracy of its kinematic data. These high-quality data 
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have well-determined orbital angular momenta and en- 
ergies, making them well-suited to the approach using 
integrals. We find patterns in the distribution of energy 
and angular momentum, never observed before, which we 
will relate to different formation and evolution processes. 
The location of our sample (its stars have a median dis- 
tance of 1 kpc) influences the processes we can study. 
The solar neighborhood lies within the inner halo 1 but 
it is also centered on the disk. We will focus in this 
paper on processes which either form stars in the inner 
halo or transfer stars into this region. We will seek to 
distinguish between stars formed very early in the inner 
galaxy (its original inhabitants) and the fashionably late 
arrivals which were accreted later from the outer halo. 
These could either have orbits sufficiently eccentric to 
bring them into the inner halo or had their orbits modi- 
fied by dynamical friction on the parent satellite to move 
them from the outer to the inner halo. 

The plan of this paper is as follows. In Section 2 we de- 
scribe our sample and investigate the properties of these 
stars in angular momentum, energy and metallicity; and 
discuss the trends that we discern. In Section 3 we will 
relate our results to the work of others, and in Section 4 
we will discuss the different formation paths for the inner 
halo, and how our data constrain them. 

2. ANALYSIS OF LOCAL HALO SAMPLE 
2.1. Sample 

Our sample is composed of red giants, red horizon- 
tal branch (RHB) stars and RR Lyrae variables with 
[Fe/H]< —1.0 and distances less than 2.5 kpc. The red 
giants and RHB st ars are all taken from the sample o f 
well-studied stars in lAnthonv-Twarog fc Twarogl (1994), 
and comprise all the stars in that study with available 
velocity data. Its RR Lyrae variables include a large 
number with new, accurate radial velocity measurements 
(van Duyne et al, in preparation) and a few with veloc- 
ities from the literature with similar or higher accuracy. 
The sample contains 246 stars, and can be accesse d at 
http: / / astronomy.case.edu/heather/fashlate.sample] 

The sample h as a large ove r lap w ith the "combined 
data sample" of iKeplev et alj (|2007f k but we have re- 
stored 12 stars which they removed from their sam- 
ple because they were likely metal-weak thick disk 
stars, in order to make the sample selection criteria 
clearer. We removed 22 giant stars which were in 
the Kepley et al sample b ut no t in the sample of 
lAnthony-Twarog fc Twarogl (|1994h . as they have dis- 
tances of significant ly lower quality. We have also added 
a few giants from lAnthonv-Twarog fc Twarogl (|1994l ) 
which were not in the Kepley et al sample because they 
had radial velocities measured after 20 00, and so were 
not in the catalog of iBeers et al.l (|2000f ) on which it was 
based. 

Our sample contains stars with accurate distances and 
well-quantified distance errors, allowing us to calculate 
errors on derived quantities like energy and angular mo- 
mentum via a Monte-Carlo procedure. There is a full 
description of our distance error calculation (which used 
a Monte-Carlo procedure to take account of all the differ- 
ent contributing errors) in Section 2.2.3 of IKeplev et al.l 

1 Contrast t he Sun's R ac of 8 kp c with distant halo stars at more 
than 100 kpc IClewlev et al.ll200a) 



(|2007f ). We emphasise the importance of our small dis- 
tance errors here: energy and angular momentum are 
strongly dependent on distance errors since two out of 
three components of velocity are obtained by multiply- 
ing the proper motion by the distance. 

Except for the three stars with new velocities, our 
sample is a subse t of the compilation of halo stars of 
IBeers et"al"1 (|2000h . We have improved the data in sev- 
eral ways. First, we removed the effect of the large 
distan ce errors in the sample found by IKeplev et al.l 
(|2007| ) by limiting the giants and RHB stars in our 
sample to those with very well-defined luminosities and 
metallicities, by using only those from th e sample of 
lAnthonv-Twarog fc Twarogl (|1994l ). The IBeers et al.l 
(|2000f ) sample contains 736 stars with [Fe/H] < — 1.0 and 
distances less than 2.5 kpc which have proper motion 
measurements. 435 of these stars are red giants, RHB 
stars or RR Lyrae variables. We estimate that 60% of 
these stars have the well-determined distances and well- 
quantified distance errors that are required for our study; 
the other 40% have data of significantly lower quality, ei- 
ther because of the lack of a luminosity measurement 
which separates stars below the horizontal branch from 
those on the horizontal branch, or because the metallicity 
measurements are less accurate and not on a consistent 
scale. However, RR Lyrae variables have large pulsa- 
tional velocity amplitudes and so require significant ef- 
fort to constra in their systemic v elocities; many of the 
RR Lyraes in Be ers etaLl (|2000l ) have large errors on 
their velocities. The restriction of our sample to those 
with both good distance and velocity measurements cuts 
the sample down to 246 stars 2 . It has a median distance 
of 1 kpc and a median distance error of 7% (error esti- 
mat es for red giant and RHB distances are given in Table 
1 of IKeplev et ^121)01 . 

Since our giant and RHB sample is selected on metal- 
licity only, and the RR Lyrae variables are selected via 
their variability, we have no kinematic selection effects 
in the sample. RR Lyraes are generally thought to be 
very old, while the giants and RHB stars could possibly 
have a larger range of age, including both younger and 
older objects. The RR Lyrae sample has a somewhat 
higher median metallicity than the giant and RHB sam- 
ple because the latter was identified with objective-prism 
searches which tend to favor extremely weak-lined stars. 

2.2. Angular Momentum plots 

iHelmi et al.1 (|1999t ) used the angular momentum com- 
ponents L z and L± to show that around 10% of local halo 
stars were formed via destruction of a satellite whose de- 
bris now occupies the inner halo. They compared the 
amount of rotation of a given star's orbit (L z ) to Lj_, 
which combines the other two components of angular 
momentum (L±— \f\L\ + Ly)) 3 - For a local sample, 
L± is dominated by v z , the velocity perpendicular to the 
plane, since y and z are close to zero, while x is around 
8 kpc. Figure [1] shows this clearly. Thus the L Z ,L± dia- 

2 We have discarded one star, HD220127, which has 
large errors on its derived velo cities using the distance from 
Anthonv-Twarog & Twarog ( 1994) and a significantly different dis- 
tance in [Schuster et al. ( 200a). 

3 L x = y v z — zv v ; L v = zv x — xv z \L z = xv y — yv x ; see the 
Appendix of Kcplcv et al. (2007) for a detailed discussion of our 
coordinate system. 
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gram is similar to the vs v z diagrams which contrast 
rotational angular momentum with the vertical extent of 
an orbit. Stars with L z less than zero are on retrograde 
orbits. 




Lperp {kpc km/s) 



Fig. 1. — Plot of angular momentum Lj_ vs the absolute value 
of the star's W velocity. 

Figure [2] shows the distribution of stars in our sam- 
ple in L z and L±. To show the location of stars with 
near-circular orbits which never reach f ar from the plane, 
the stars from the iNordstrom et alj (|2004h sample, a 
magnitude-limited sample of stars with distances up to 
^200 pc, are shown in grey on the same diagram. The 
Nordstrom et al. sample is almost completely dominated 
by thin disk stars. As expected for a sample of stars in 
the thin disk whose orbits are confined to the plane, few 
of these stars have L± greater than 200 kpc km s _1 , and 
their L z values are higher than almost all stars in our 
h alo sample. 

iKepley et al.l (|2007h have di scussed star s tream s in the 
sample: outliers such as the iHelmi et al.l (|1999| ) group 
near (X Z ,L^)=(1000,2000) kpc km s _1 and the very ret- 
rograde stars with (L Z ,L±) near (-2000,1500) kpc km 
s _1 . In this paper we focus, not on these outliers, but 
on the majority of the stars in the sample which occupy 
the region normally described as the smooth, well-mixed 
halo. Thus, for the rest of the paper, we consider stars 
with L±< 1400 kpc km s _1 . Figure [3] shows a contour 
plot of the stellar density in this diagram. 

The classical description of halo ki nematics uses a ve- 
locity ellipsoid (jSchwarzschildl 119071 ) with independent 
Gaussian distributions in v x ,v y and v z . In the small 
range of distance in our sample (recall that the median 
distance is 1 kpc) these assumptions will lead to roughly 
Gaussian distributions of L z and L± as well. It is clear 
from Figure [2] that such a velocity ellipsoid is not a good 
description of our sample, even when outliers are re- 
moved by restricting it to stars with L±< 1400 kpc km 
s _1 . The only part of the sample that appears smooth 
and close to Gaussian is the region occupied by stars 
whose orbits are quite confined to the plane: those with 
low L± (< 400 kpc km s _1 ). This is the region of L± 
occupied by disk and thick disk stars in Figure[2] For val- 
ues of L± above 400 kpc km s _1 , the distribution looks 
quite clumpy. For example, there is a pronounced gap 
around (L Z ,L±) = (500,400) kpc km s^ 1 and clumps at 



(750,500) and (0,900) kpc km s _1 . 

2.3. Statistical Analysis 

We use two statistical approaches to quantify the struc- 
ture in the [L Z ,L±) diagram. In the first approach, which 
we present here in detail, we fit mixtures of Gaussian 
distributions to L z distributions for different ranges of 
L± , decide how many components are needed in each L ± 
range (based on statistical model selection procedures), 
and then estimate the value of L± where the distribu- 
tions change character. Then we compare the results 
of the first mixture fits with those based on the actual 
estimated "change points" (the values of L± where the 
distributions change character). 

We begin by dividing the sample by eye into three 
ranges of L±: (0-350,350-700,700-1400) and determin- 
ing whether each distribution of L z is best described by 
a single Gaussian: 

or a mixture of two: 

/ ~ pNfaual) + (1 - p)N(fi 2 , a 2 2 ), < p < 1 

Figure 2] shows the three histograms. It can be seen 
that the histogram for the lowest L± values (stars whose 
orbits are confined to the disk plane) looks unimodal, 
while the other two (stars with orbits which reach further 
from the plane) appear less so. 

We decided whether to use a unimodal or bimodal 
model using an extension to the standard procedure of 
calculating the likelihood of the sample under a given 
model, the Akaike Information Criterion (|Akaikelll97l 
AIC hereafter): 

AIC = -21ogL(0) + 2d 

Here L{9) is the likelihood function for the model with 
parameters 9 and d is the number of parameters in a 
competing model, unimodal or bimodal. A model with 
one extra parameter is thus penalized by the factor of d 
when comparing values of log(likelihood). 

Model comparisons using the AIC can be seen in Ta- 
ble [H where the preferred models are shown in bold 
face. A competitive model selection criter ion is Schwartz ' 
Bayesian Information Criterion (BIC, Schwarz Il978f) 
which changes 2d to log(n).d, where n is the number of 
data points. If we had used the BIC for the L±= 350 
range, it would have favored the one-component model 
even more strongly, because this r ange of L± conta i ns the 
largest number of stars. However. Ilshwaran et al.l (|2001l ) 
have demonstrated that AIC is a better criterion for com- 
parison of mixture models, so we prefer to use this for 
our analysis. It should be noted that although the differ- 
ence between the AIC values for one and two components 
looks small compared to its large value, a difference of 2 
in its numerical value is equivalent to the amount which 
justifies an additional parameter; also, different analyses 
discussed below give similar results. 

We have shown that when we divide the sample into 
three L± ranges by eye, we find that the stars whose 
orbits arc confined near the plane have a different dis- 
tribution than those whose orbits reach further from the 
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Fig. 2. — Plot of angular momenta L z (measuring rotation) vs L± (in our sample, this correlates with the distance a star's orbit reaches 
from the disk plane) for stars with [Fc/H]< — 1. It can be seen that the er rors on these quantities are quite small for most of the sample. 
RHB stars are plotted with open circles. Stars from Nordstrom et al. (2004) with [Fc/H]>-1.0 (which are predominantly from the disk and 
thick disk) are plotted in grey, and the globular cluster M4 is shown with an open 5-pointed star. 
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Fig. 3. — Kernel Density Estimate plot of angular momenta L z 
vs L± for stars in our sample with [Fe/H] < —1.0. 



TABLE 1 

Model Comparison for data grouped in 

L± BY EYE 



Lj_ range 


n 


AIC(l) 


AIC(2) 


0-350 


103 


1685.431 


1686.274 


350-700 


63 


1041.645 


1035.361 


700-1400 


61 


1018.468 


1015.161 



plane: the former have a simple Gaussian distribution 
of L z , while the latter have a bimodal distribution. We 
can now remove the effect of our subjective choice of L± 
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Fig. 4. — Histograms of L z for stars with [Fc/H] < —1.0, for 
three ranges of Lj_ as shown. 

range using a change point analysis (ICsorgd and Horvathl 
H997L [Chen and Guptal l2?M lGanocvll200.^ . which will 
solve for the best values of L± to separate the three 
groups. In particular we are interested in the value of 
Lj_ where the distribution changes from unimodal to bi- 
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modal, and in interpreting this number in terms of the 
Galaxy's known stellar populations. 

For the change point analysis we fit the following model 
to the data, solving for the two dividing points and 

fi(L z \L ± <L ±il )~NQi u al) 
h(L g \L ±il <L ± < L ± , 2 ) ~ P N(fi 2 ,a 2 2 )+{l-p)N{fi 3 ,a 2 3 ) 

f S (L z \L ±>2 <L ± < 1400) ~ N(^aj) 

, where < p < 1. 

(We have chosen a single Gaussian for the highest L± 
group to make the change point analysis simpler; note 
that the middle group has the strongest deviation from 
Gaussian shape.) 

The dividing points which produce the highest like- 
lihood are (L± t %, £±,2) = (349,709) kpc km s _1 ; these 
are quite close to the original values chosen by eye. Re- 
peating the analysis summarized in Table Q] with these 
slightly different values produces the same result: the so- 
phisticated AIC criterion prefers a Gaussian distribution 
for L± less than 349 kpc km s _1 , and mixtures of two 
Gaussians for the other two ranges of L±. 

The lower change point, L±= 349 kpc km s^ 1 , is par- 
ticularly interesting because it represents the value of L± 
below which the distribution of L z is Gaussian. A rough 
calculation (setting y = z = and x = 8) shows that this 
L± "change point" corresponds to \v z \ ~ 40 km s^ 1 . The 
z velocity dispersion of all stars in our sample with L±< 
350 kpc km s _1 is 30 km s -1 , 4 intermediate between the 
thin disk (a z = 20 km s _1 ) and the thick disk (a z — 40 
km s _1 ). 

Table [2] shows the best fit mean and sigma for these 
models. First, it should be noted that none of the groups 
have rotational properties at all close to the thin and 
thick disks' «£ z }~ 1700 kpc km s" 1 ). The low L± 
group (stars with orbits closely confined to the plane) 
shows a mild prograde rotation ((L z )= 364 kpc km s _1 , 
corresponding roughly to = 45 km s _1 ). The other 
two groups (containing stars with orbits reaching further 
from the plane) have more complex distributions of Lz, as 
can be seen in Figure 2J Lz is roughly zero in the mean 
but quite clumpy, making the mean a less informative 
parameter. 

It is important to note here that the somewhat higher 
mean rotation of the group with orbits very close to the 
plane is not due to thick disk contamination: Figure [2] 
shows that there are almost no stars in this lowest L± bin 
in the region occupied by disk and thick disk stars, which 
have near circular orbits. The low Lj_ stars in our sample 
are on highly eccentric orbits close to the pl ane, as was 
first noted bv lSommer- Larsen fc Zhenl (|1990l) . The pres- 
ence of this low scale height component of the metal poor 
halo was not found by previous studies, perhaps because 
of the significantly larger errors or the high-latitude sam- 
pling of some surveys which led to a lack of data close to 
the plane. 

4 Fortuitously, these two quantities do not vary exactly in lock- 
step, which would mean that truncating Lj_ would produce a tru- 
cated distribution in v z and make estimation of the velocity dis- 
persion difficult. 



What is the spatial distribution of these stars with low 
L±7 We cannot measure it directly from our sample, be- 
cause there were latitude-related selection effects in the 
initial surveys on which it was based (see Section 3.2 be- 
low), and our sample contains only stars from these sur- 
veys which have been followed up extremely thoroughly 
by a number of authors, thus giving very accurate kine- 
matical measurements, rather than all stars identified 
down to a certain magnitude by the surveys. It is also not 
possible to simply compare the z velocity dispersion of 
these stars with that of the thick disk, because the thick 
disk's large degree of rotational support adds extra flat- 
tening. However, we can make an estimate of the axial 
ratio of this component using its v elocity ellipsoid, and 
emplo ying the tensor virial theorem. iBinnev &: Tremaind 
(1987) discuss the case of a self-gravitating population 
in the contex t of the flatte ning and rotation of elliptical 
galaxies, andlHclmi (2008) discusses the case of a tracer 
population such as the stellar halo. We find that the ve- 
locity ellipsoid for stars with L± less than 350 kpc km 
s _1 , {<jjj , ay , aw) = (144,111,30), corresponds to a flat- 
tening c/a = 0.2. This contrasts markedly with the flat- 
tening for the entire halo at the solar radius (e/a=0.6), 
which can ei ther be derived direct ly from the models of, 
for example, iPreston et alj (119911) or from the local ve- 
locity ellipsoid (see Hclmi 2008). 

In our second statistical analysis, we used the model- 
based clustering method in the R project for statistical 
computing 5 package MCLUST 6 . In the MCLUST anal- 
ysis, the data are fit by a mixture of a finite number 
of two-dimensional Gaussian distributions, and the best 
fit is computed using a modification of the AIC with a 
second order correction for small sample sizes (AICc 7 ). 
We considered both 2-D Gaussians with major and mi- 
nor axes parallel to the L z and L± axes, and more com- 
plex models where the Gaussians could be oriented in 
any direction. A good fit was obtained with two ellip- 
tical Gaussians. The one at low Lj_ is quite flattened, 
with major axis parallel to the L z axis; the other, for 
the higher Lj_ values, is a rounder elliptical Gaussian. 
These Gaussians can be seen in Figure \5\ The separa- 
tion point between the two elliptical Gaussians is close 
to L z = 350, confirming, with a different analysis tech- 
nique, the results of the change point analysis reported 
above. Another analysis, using the AIC only, divided the 
data with L± greater than 350 into five different clusters, 
so we should not take the two-Gaussian result shown in 
Figure [5] to indicate a lack of substructure in the data 
with high L±. 

2.4. Binding Energy 

It will be useful in our discussion of the different pro- 
cesses which deposit stars in the inner halo to consider 
the total energy of each star, as well as its angular mo- 
mentum. Figure [6] shows the traditional energy-angular 
momentum plot. (Note that we have computed these 
values for the entire sample, not just those stars with L± 
less than 1400 kpc km s _1 .) To ca lculate total energy we 
have used the Galaxy potential of iJohnston et all (|1995l ) 
with a somewhat smaller disk mass of 7.5 x 10 10 M©and 

5 http://www.r-project.org 

6 http:/ /cran.r-project.org/doc/packages/mclust.pdf 

7 http:/ /en. wikipedia.org/wiki/Akaike_information_criterion 
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TABLE 2 

Parameters of best fit Gaussian mixtures 
for l z distributions for different l± 

RANGES 



L ± range p fii a\ ^2 &2 

0-350 1.0 364 848 

350-700 0.45 -460 348 1088 567 

700-1400 0.70 186 1035 54 194 




-2000 -1000 1000 2000 3000 

Lz (kpc km/s) 



Fig. 5. — MCLUST fit of elliptical Gaussians to our data. It 
can be seen that this alternative analysis produces very similar 
answers: a highly flattened elliptical Gaussian for the low values 
of Lj_, separated from the other data near Lj_= 350 kpc km s _1 . 

a dark halo circular velocity of 155 km s _1 . We set the 
potential energy equal to zero at R=200 kpc. 

It can be seen that the majority of stars in our sam- 
ple have higher binding energy (lower Etot) than the en- 
ergy corresponding to the LSR velocity. This is expected 
in a sample of halo stars because of the halo's steep 
density distribution (i? -35 , measured for RR Lyraes by 
I Vivas fc Zinnl l2006). The dumpiness that we see in the 
L Z ,L± diagram is also visible in this diagram, particu- 
larly along the lines L z =1200 and kpc km s _1 . The 
clump near L 2 =1100 kpc km s _ \.Et ot =-160,000 km 2 /s 2 
has energy a nd angular momentum similar to the Arc- 
turus group (|Eggenlll998l iNavarro et al.ll2004f ). We also 
see the "plume" at L z ^ —300 kpc km s~ 1 an d a large 
range of E to t which has been remarked upon by iDinescvJ 
( 2002) and lBrook et~aT1 (|2003h . They suggested that this 
"plume" is debris from the progenitor of u> Centauri. 

In Figures [7] and [8] we see the intriguing result that the 
RR Lyraes and red giants in our sample show differences 
in binding energy 8 . (Note that we have excluded the 
RHB stars from this comparison for clarity.) The red 
giants (which have significantly more metal-poor stars) 
are on average more tightly bound to the Galaxy than the 
RR Lyraes. We find that 55% of the RR Lyrae sample 
have binding energy less than the LSR's, compared with 
only one third of the red giants. We tested the statistical 
significance of this result by sampling the (122 star) red 
giant sample with replacement 1000 times in order to 
produce samples of the same size as the RR Lyrae sample 
(99 stars) . None of the thousand samples had more than 
55% with binding energy less than the LSR's, and only 

8 The differences between the two samples in the L z -L^ diagram 
are much more subtle. 
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Fig. 6. — The energy-angular momentum diagram for our sam- 
ple. Points are coded by m etallicity as in Figure [2] and the pre- 
dominantly disk stars from Nordstrom et al. (2004) are shown in 
grey. RHB stars are again plotted with open circles. Black lines 
with closed circles show the circular orbits at a given position, with 
orbits at 6 and 8 kpc marked, while the parabola at the bottom 
shows the position of orbits in the plane with apogalacticon 7 kpc. 
The star is the globular cluster M4, which has a near-radial orbit 
and is currently 6.3 kpc from the galactic center. 

one had half the sample with less binding energy than 
the LSR. It is clear that the giants and the RR Lyrae 
variables have different distributions in binding energy. 
We note that since less than 10% of our sample are RHB 
stars, it is not possible to make a similar check for these 
stars. 

This result is a kinematic counterpart of the 
iPreston etail (|1991[ ) result that RR Lyraes and BHB 
stars have different spatial distributions, with the RR 
Lyraes less centrally concentrated than the BHB stars. 
It is also related to the different population properties 
of both field RR Lyraes and globular clusters when they 
are divided according to the Oosterh off dichotomy (eg 
iLee fc Carneylll999llMiceli et al.|[2008l) which is likely re- 
lated to age (|Lee fc Carnevlll999l . [Jurcsik et al.ll2003[ ). 

The giants in our sample will evolve onto the hor- 
izontal branch as they age, and since the horizontal 
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Fig. 7. — The energy-angular momentum diagram for our sample, showing RR Lyrae variables on the left and (first ascent) red giants 
on the right. Points are coded by metallicity as in FigurefS] Lin es mark loci of circul ar orbits and apogalacticon 7 kpc, as in Figure[6] The 
orange star is the globular cluster M4, which has [Fe/H]=-1.22 (Kraft & Ivans 2003). 

branch morphology in the solar neighborhood is very blue 
(|Kinman et al.|[l994f ). most of them will become BHB 
stars. Thus we see that the lower binding energy of our 



RR Ly raes is consistent with the result of iPreston et al.l 
(|1991h that RR Lyraes have a shallower spatial distribu- 
tion than BHB stars, and so are more likely to be found 
in the outer halo. Our RR Lyrae sample contains more 
objects from the outer halo than the red giant sample, 
although not all RR Lyraes in our sample arc outer halo 
objects, as can be seen in Figure [H 




-15 -10 
E lol (10 4 kmVs 2 ) 



-2000 2000 
Lz (kpc km/s) 
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Fig. 8. — Histograms of Etot L z and [Fe/H] for the RR Lyrae 
and red giant subsamples. RHB stars have been omitted. The 
energy of the LSR orbit is shown on the left panels with a dotted 
line. It can be seen that the giants are more tightly bound on 
average than the RR Lyraes. 



Connecting horizontal-branch morphology to age 



would suggest that these outer h alo objects we r e also 
younger, as originally su ggested bv lSearle fc Zinnl (|1978l ) 
and then by I Zinnl (jT993). The younger age of an accreted 
outer halo compared with the inner halo formed in situ 
makes sense in the light of the theoretical expectation 
that the first stars to form in the Galaxy would be in the 
rare hi gh-density peaks, and would be centrally concen- 
trated (|Diemand~etI^ EoqII) ; outer halo objects would 
form somewhat later. The dominance of retrograde or- 
bits in the low binding-energy stars of bot h types is also 
intere sting in the light of the discovery of ICarollo et al.l 
(|2007f ) that the outer halo has retrograde rotation. It 
will be very interesting to see whether the nearby BHB 
stars (which may be older than the RR Lyraes) share the 
properties of the red giants in our sample in this diagram. 

There are other interesting differences between the RR 
Lyrae and red giant kinematics in our sample: many 
of the more bound red giants are have prograde orbits, 
including a number close to the circular orbit line, while 
the RR Lyraes have a smoother distribution of orbital 
angular momenta. 

2.5. Trends with Metallicity 

Trends with metallicity are interesting because they 
give a handle on two processes: dissipational collapse 
with star formation subsequently enriching later gen- 
erations, and dynamical friction, which is more effec- 
tive for massive objects which produce stars with higher 
mean m etallicity becaus e of the mass-metallicity rela- 
tion (eg iLee et all |2006|) . We see an interesting differ- 
ence in kinematics between the very metal-poor stars 
(with [Fe/H] < —1.5) and the metal-richer halo stars 
(with —1.5 < [Fe/H] < -1.0), in Figure ® where stars 
with different metallicity are shown in different colors. 
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While the entire sample has values of L z ranging from 
-3000 to 2700 kpc km s _1 , there are no stars with in- 
termediate metallicity (— f.O > [Fe/H] > —1.5) which 
have L z below -800 kpc km s _1 . (Note that with a m e- 
dian halo metallicity of -1.6, eg lRvan fc Norriill99lL , 
these intermediate-metallicity stars account for almost 
half of the halo, although they are under-represented in 
our sample because of the relative ease of identifying very 
metal-poor giant stars in surveys.) 

We see the trend with metallicity more easily in the 
histograms of Figure [9l The median L z value for the 
metal-richer stars is 545 kpc km s _1 , significantly higher 
than the median fo r the metal-poor stars (135 kp c km 
s _1 ). The sample of iSommer-Larsen fc Zhenl (1990) also 
shows this difference in L 7 for more metal-rich halo stars. 




L x (kpc km/s) 

Fig. 9. — Histograms of L z values for stars with L±_ less than 
1400 kpc km s~ x . Upper panel shows stars with [Fe/H] < —1.5, 
lower panel those with [Fe/H] between —1.0 and —1.5. 



3. DISCUSSION 

3.1. Halo flattening: another component 

The first suggestions that the halo's flattening might 
be variable wi th r adius were in the 19 50s and 60s by 
ISchmidtl (| 1956( 1 and lKinman et alJ (j!965h . studying glob- 
ular clusters and RR Lyraes respectively. While both 
studies were contaminated by some metal-richer objects, 
it is u nlikely th a t this was the only reason for their find- 
ings. lHartwickl (|1987l ) suggest ed a t wo-co mponent halo 
based on RR Lyrae samples. IZinnl (|1993f) also found a 
two-component halo by studying the kinematics, space 
distribution, metallicity and ages of globular clusters in 
the Galaxy. The Old Halo clusters of the inner halo show 
a flattened distribution and have some rotation (although 
much less than the disk's). The Younger H alo is spheri- 
cal an d has no net rotation. More recently, iMiceli et al.l 
( 2008) discussed a two-component halo with distinct ra- 
dial distributions for different Oosterhoof types in their 
large sample of RR Lyraes, and lCarollo et all (|2007f ) used 
a very large sample of halo stars from the SDSS survey 
to extend the description of the two-component halo. Of 



particular importance to an understanding of the popu- 
lation characteristics of the two components is Carollo et 
al's demonstration that the near-spherical outer halo has 
a net retrograde rotation and that it is more metal-poor 
than the flattened inner halo. 

All of the field star results are for relatively distant 
samples. What of samples like o ur own, with many stars 
less th an 1 kpc from the plane? ISommer-Larsen fc Zhenl 
(1990) used an early local sample of 151 metal-poor, pre- 
dominantly giant and RR Lyrae stars with less accurate 
data than ours. Using only stars with [Fe/H] < -1.5, 
they reconstructed the density distribution of the halo, 
found that these metal-poor stars had a two-component 
structure, and noted that their flattened shape was due 
to velocity anisotropy, not rotation. While they do not 
explicitly quote the axial ratio of the flattened compo- 
nent, it can be seen in their Figure 6 that it is quite flat: 
c/a around 0.25. In general, our analysis agrees wel l with 
these early results. ISommer-Larsen fc Zh cn ( 1990;) sug- 
gested that either the metal-poor stars in the flattened 
component formed in an anisotropic, dissipative collapse 
(with the anisotropic infall caused by a flattened dark 
halo) or by late infall of metal-poor gas clouds which 
" plunged into the disk " and formed stars there. 

iPreston et al.l (|1991f ) used BHB stars and RR Lyraes 
(which were dominated by stars more than 2 kpc from the 
plane except for some very close to the galactic center) 
and derived a density law for the halo which became more 
spherical with increasing radius, with axial ratio c/a = 
0.5 at the galactic center. However, when they used this 
model to predict the RR Lyrae density near the Sun, 
they found that it under-predicted the known value by a 
factor of two, concluding that "there exists a metal-poor 
population, confined to a volume near the galactic plane, 
which has remained undetected in intermediate and high- 
latitude surveys." It seems very likely that Preston et al 
were descr ibing what we call the lo w L± component of 
the halo. I Green fe Morrison] (|1993l l also suggested the 
existence of this highly flattened halo component using 
a sample of BHB stars within 500 pc of the Sun. 

Thus there are two flattened components in the in- 
ner halo: the m oderately flattened, c/a ~ 0.6 compo - 
nent identified bylKinman et all dl 965ft jHartwickl (11987ft 
IPreston et all (l 99 1( 1. IZinnl (|1993l ) and iKinman et all 
(| 19941 ) from more distant samples of halo stars, and 
th e highly flattened, c/a ~ 0.2 component foreshadowed 



i com pc 

bv ISommer-Larsen fc Zhenl Jl990( ). IPreston et all (|1991l ) 
and iGreen fc Morr ison (1993) and revealed in more de- 
tail in this work. We will discuss possible origins for these 
t wo components below. 

IChiba fc Yoshiil (fl998ft and lChiba fc Beersl (f2000h also 
studied local samples, although the one used by the lat- 
ter (the sample of Beers et al. 2000) suffered from large 
distance errors for a number of stars (see S ection 2.1 and 
iKeplev et al.1l2007f l. IChiba fc Yoshiil |l998ft focused their 
attention on quantifying the contribution of the metal 
weak thick disk to t heir sample. Curiously, although 
IChiba fc Beersl (|2000) use a local sample, they conclude 
that their results are similar to those of the more dis- 
tant samples such as IPreston et al.l (|1991[ ): a halo with 
an axial ratio decreasing as R decreases and only mod- 
erate flattening in the center. We suggest that the large 
distance errors of some of the stars in their sample might 
have made it difficult to distinguish between a moder- 
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ately flattened and a highly flattened inner halo. We 
explored the effects of a larger measurement error us- 
ing our data, by adding a random number picked from a 
Gaussian with twice the error in L z and L± to each dat- 
apoint. Much of the structure that we see is no longer 
visible, showing that our careful sample selection was 
needed. 

Recently, I Brown et al.l (|2008f ) have presented a prelim- 
inary kinematic analysis of a large sample of 2414 BHB 
stars with distances larger than 2 kpc from the galac- 
tic plane. They find a flattened component, but identify 
it with the metal-weak thick disk because of its signifi- 
cant rotation. However, in their paper they simply solve 
for mean rotation of the entire sample as a function of 
distance from the plane, and find both a trend with z 
and a surprisingly large random variation in mean ro- 
tation (from +150 to -60 km s" 1 ) over the ran g e of z 
heights covered by their sample. iKinman et alj (2008) 
have shown, using proper motion information for the 
most nearby stars in this sample, plus another sample 
of nearby BHB stars, that the mean rotation of BHB 
stars close to the plane i s close to zero Th us the flat- 
tened component in the iBrown et al.l (|2008h sample is 
likely to be associated with our low L± group. 

It is also interesting to note that the nearest globu- 
lar cluster, M4, has properties consistent with member- 
sh ip in the highly flatt ened halo. Using the data given 
in iDinescu et alj (|1999f ) we find that it has z=-0.b kpc, 
L z = 215 ± 141 kpc km s _1 and L x = 32 ± 22 kpc km 
s _1 . We show its position in the energy and angular 
momentum diagrams of Figures [5] and El with an open 5 - 
poi nted star. Using the data fro m IDinescu et al.1 (|1999f ) 
and lCasetti-Dinescu et al.1 (|2007f l. we find that there are 
another four globular clusters with R gc between 6 and 10 
kpc which have similar orbital properties to our low L± 
group: NGC 4372, 4833, 5139 (u Centauri) and 6779, 
and another four with R gc less than 6 kpc: NGC 5986, 
6093, 6656 and 6712. 

3.2. The thick disk 

It is remarkable how little contribution the thick disk 
makes to our sa mple; the region o f the L z -L± diagram 
occupied by the iNordstrom et al.1 ((2004) local dwarfs, 
a sample almost totally dominated by the disk, con- 
tains less than 5% of our sample. Given the fact that 
thick disk stars outnumber halo stars by 1-2 orders 
of magnitude in our region, this could put strong con- 
straints on the size of the low-metallicity tail of the thick 
disk metallicity distribution . Earlier claims of a sig- 
nificant low-metallicity tail (|Morrison etal]ll990h were 
over-estimates becau se of problems with the DDO metal- 
licity calibration (iTwarog &: Anthony- Twarod I1994I . 
iRvan fc Lambertlll995l IChiba fc Yoshiilll998f ). 

However, both the RR Lyrae and the red giant/RHB 
samples are seriously incomplete at low latitudes, as can 
be seen in Figure [TU1 This incompleteness would reduce 
the number of thick disk stars found in our sample. While 
this is unlikely to be a strong effect for our very local sam- 
ple (recall that the mean distance of our stars is 1 kpc), 
we cannot use this sample to obtain an unbiased estimate 
of the density of metal-weak thick disk stars unless we 
engage in complex correctio ns for our se le ction effects. 
By contrast, the sample of iNorris et alj (|1985l ) which 
was used for first studies of the metal-weak thick disk 




lamude [b| 



Fig. 10. — Histograms of latitude \b\ for (top) the RR Lyraes in 
our sample (middle) the r ed giants and RHB stars in our sample, 
and (bottom) the giants of Norris et al. (1985j). The curve on each 
plot is cos(b). It can be seen by comparison with the curve on 
each plot that our local sample has serious incompleteness at low 
latitude, and the NBP sample is much more complete. 

(|Norris et al.| [T985. Mor rison et al .119901 ) shows much less 
incompleteness with latitude, as can be seen in the bot- 
tom panel of Figure 6. An upcoming analysis of the 
metallicity of t he stars in this sam ple (Beers, private 
communication, iKrugler et aTll2004D would be the best 
way to obtain an unbiased estimate of the local density 
of the metal-weak thick disk. 

3.3. Summary of Results 

Pulling together all of the above analysis, we have the 
following results. 

(1) We find evidence for an additional component to 
the inner halo, which is highly flattened (c/a ~ 0.2), 
pressure-supported and has a small prograde rotation. 
This component is not associated with the metal-weak 
thick disk, which is predominantly rotationally sup- 
ported. The new component shows a smooth distribution 
of L z . 

This component is separate from the moderately flat- 
tened inner halo which was discovered using samples of 
halo stars l ocated further t han 1 kpc from the Sun , for 
example by IKinman et ail (|1965l ) , lHartwickl (|1987f ) and 
iPreston et al.l (|1991[ ). This moderately flattened compo- 
nent has axial ratio (c/a) of ~0.6. If the globular clusters 
are representative of the field stars in this moderately 
flattened component, it is pre dominantly o l d and also 
has a small prograde rotation. ICarollo et al.l (2007) also 
ascribe a "modest prograde rotation" to their moderately 
flattened inner halo. 

The highly flattened component forms about 40% of 
our sample. However, this could be an under-estimate 
of its true density as the samples on which we base our 
analysis avoid low latitudes (see Figure IT0|> . 

(2) The remaining 60% of stars in our sample, whose 
orbits reach further from the plane, have a mean rota- 
tion close to zero and a clumpy distribution in angular 
momentum and energy. 

(3) We see trends with metallicity. The low metallicity 
stars in our sample ([Fe/H] < —1.5) show a mean rota- 
tion near zero and a clumpy distribution in energy and 
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angular momentum. 

The intermediate-metallicity stars (—1.5 < [Fe/H] < 
— 1.0) cover the whole range of L± (in other words, they 
belong both to the highly flattened component and to 
the other (s)), but there is a higher proportion of these 
intermediate-metallicity stars in the low L± region. The 
intermediate-metallicity stars have a small prograde ro- 
tation and a clumpy distribution in angular momentum 
and energy. 

(4) The RR Lyrae variables in our sample have a lower 
binding energy, on average, than the red giants, with half 
of the RR Lyraes being less tightly bound than the LSR 
orbit, in contrast to the red giants where only one third 
of the stars are less tightly bound. Because the local 
horizontal branch morphology is very blue, this suggests 
that the younger RR Lyrae stars are less strongly bound 
to the Galaxy. The RR Lyrae sample also appears to 
contain more stars on retrograde orbits. 

4. FORMATION AND EVOLUTION OF THE INNER HALO 

We will now discuss the different stages of halo forma- 
tion and the constraints placed on various theories by our 
data. In particular we will examine possible formation 
paths for the highly flattened halo component which is 
supported by velocity anisotropy. 

4.1. The First Halo Stars 

Theory predicts that the first stars to form in the 
Galaxy's halo occupy dark halos which formed from very 
rare high-density peaks in the primordial density field. 
At this early stage (perhaps beginning earlier than z=10, 
less than 0.5 Gyr after the Big Bang) t he Galaxy's disk 
had not formed. iDiemand et al.l (|2005f ) use simulations 
to predict that the density distribution of such stars is 
more centrally concentrated than the rest of the halo. 
Many of these first stars will now be found in the inner 
few kpc of the Galaxy; however there will still be a few 
with orbits that reach the solar neighborhood. 

How might we identify these first stars? We would 
expect the very first stars to have extremely low metal- 
licity. However, the amount of chemical enrichment that 
occurs in the earliest star formation in the halo is not well 
constrained, so we do not know whether only extremely 
metal-poor stars were formed at this stage, or whether 
some stars closer to typical halo metallicities were also 
formed. We see no unique pattern in the angular momen- 
tum or energy distribution of the most metal-poor stars 
in our sample. We plan to extend our sample to stars of 
extremely low metal abundance and check whether their 
kinematics are different. 

Where would the early halo stars, formed long before 
the disk, appear in the L z -L± diagram? We would ex- 
pect them to have low total energy (be tightly bound to 
the Galaxy) and have both low L and low L z . Interest- 
ingly, the very metal-poor st ar CD-38 245 ([Fe/H1=-4.0) 
is in t he original sample of lAnthonv-Twarog fc Twarod 
(fl99l but was excluded from our sample because of its 
distance (3.3 kpc). However, its angular momentum and 
energy are well-determined, and fit the expectation for 
a star formed very early: L z — -110 kpc km s" 1 , Lj_ = 
254 kpc km s _1 and E to t= -1.6 xlO 5 km 2 s~ 2 : it is on a 
near-radial orbit, confined to the inner halo. It is quite 
possible that CD-38° 245 was formed early, and is a gen- 
uine "original inhabitant" of the inner halo. 



Because the Galaxy itself is quite small at such early 
times, the first halo stars to form would have low total 
angular momentum, and thus low L±. However, they 
cannot form our entire low Lj_ group, because of its 
large range in L z . Also, it is unlikely that the chaotic 
conditions under which the first stars formed would al- 
low a very flattened configuration, such as this group, 
to survive. This is also the reason why any disk which 
formed during the early period where the Galaxy's halo 
was growing via violent merging would not survive. 

We should note, however, that there is another evolu- 
tionary path which will also produce stars with low angu- 
lar momentum and energy: dynamical friction dragging 
satellites in from the outer halo. So the possession of low 
angular momentum and low energy are not sufficient to 
identify early halo stars. 

4.2. Accretion and halo growth before disk formation 
4.2.1. Substructure in E and L 

Before the disk forms, mergers will cause the Galaxy's 
gravitational potential to vary strongly. This violent re- 
laxation might smooth out any dynamical record of the 
origins of stars formed then, and lead to a kinematically 
smooth distribution. This has been the source of a gen- 
eral expectation that the inner halo be smoother than the 
outer halo. However, in our inner halo sample, the only 
stars which display a smooth L z distribution are ones 
which occupy the highly flattened component. We have 
argued above that we would not expect this component 
to be formed before the disk because such a flattened spa- 
tial distribution would be destroyed by mergers. Thus its 
smoothness cannot be ascribed to violent relaxation and 
is not relevant here. 

However, we would expect many of the stars in the 
remainder of the sample to be formed at this stage of 
the halo's growth, and we see no evidence for a complete 
smoothing of kinematical quantities by violent relaxation 
there. We are not aware of any other observational study 
which finds such a clumpy distribution in energy and an- 
gular momentum. We suggest that the larger distance er- 
rors of previous samples (which become magnified when 
quantities such as energy and angular momentum are 
measured) have obscured the intrinsic dumpiness of the 
data there. We note that our clumpy distribution in 
angular moment u m an d energy does not contradict the 
result of iGouldl (|2003h that no more than 5% of the 
local halo comes from a given stream because multiple 
streams can come from a single progenitor. Results from 
simulations of the dark h alo (|Helmi et al.l l2002t 120031 . 
IDiemand et al.l 2005: 2007) have a similar result: the ve- 
locity distribution of the dark matter is quite smooth 
there, despite a small number of progenitors supplying 
the field stars of the inner halo. 

The clumpy distribution in energy and angular mo- 
mentum implies that either the violent relaxation was 
not violent enough to destroy structure in angular mo- 
mentum and energy, or this process was limited to the 
inner few kpc of the Galaxy, and did not reach the solar 
neighborhood. This is good news for galactic archaeol- 
ogy, as it will allow us to trace various halo substructures 
to very early times. If, in addition to the preservation 
of some initial conditions through the violent relaxation 
stage, a small number of accreted objects formed most of 
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the star s in the inne r halo , as suggested by the s imula- 
tions of lHelmi et all (|2002fl and lHelmi eFatl (|2003l ). this 
would give a natural explanation for the dumpiness we 
see. It is unlikely that such a clumpy distribution in an- 
gular momentum would be produced by a large number 
of progenitors with distinct orbital properties. 

4.3. Pure accretion 

The corollary of the expectation that the earliest stars 
will form in the high-density peaks near the center of 
the Galaxy is that star formation in the lower-density 
regions that join the Galaxy as its halo grows should 
start somewhat later. This leads to the expectation that 
RR Lyrae samples, which should be younger on average 
than BHB samples, should be less strongly bound to the 
Milky Way. We have discussed the fact that we see this 
in our sample in Section 2.4. 

Now let us focus on "pure" accretion, unassisted by 
dynamical friction, which builds up the outer halo, and 
ask whether any of these outer halo objects would be 
found in our local sample. As the Galaxy grows, and 
its turnaround radius (where objects separate from the 
Hubble flow and become bound to the Galaxy) increases, 
it becomes harder for accreted objects to reach the so- 
lar neighborhood unless they are on near-radial orbits. 
However, radial o r bits a re quite common for inf ailing 
satellites: iBensonl ()2005f ) shows that more than half of 
the orbits of satellites in the simulations he studies have 
eccentricity between 0.9 and 1.0 when they enter their 
host halo, showing that this is not an impossible require- 
ment 9 . 

Such objects, with near-radial orbits, have all com- 
ponents of angular momentum near zero. If dynamical 
friction has not occurred, then we have the additional 
diagnostic that the orbits will have low binding energy 
(high total energy) as well, since they fell into the Galaxy 
from a large distance. Stars from these late-accreted 
small objects will be found in the region near L z ~ and 
E to t greater than, say, -1.4xl0 5 km 2 sec -2 . Intriguingly, 
there is a preference for slightly retrograde orbits in this 
region of the L z -E to t diagram, and we see more of this 
component in the (possibly younger) RR Lyrae sample 
th an in the red giants. This may be related to the result 
of lCarollo et al.l (|2007l ) that the outer halo has a small 
net retrograde motion; if the material which was added 
to the Milky Way at later times had a different mean 
angular momentum from the material which formed the 
disk, then the only part of this outer halo that could 
reach the solar neighborhood without help from dynami- 
cal friction would be the stars on near-radial orbits. The 
remainder of these stars would then populate the outer 
halo. However, it must be remembered that dynamical 
friction can also populate this area, and in fact many 
of these stars have als o been claimed to be the debris 
of uj Cen's proge nitor (|Dinescull200l iBrook et al.ll2003l 
iMeza et al.ll200l . 

We have noted that it is not likely that our highly 
flattened halo component formed before the disk be- 
cause the ongoing mergers which only stopped when the 
present-day disk formed would tend to heat this very 
flattened component. However, the moderately flattened 



inner halo could have been formed during mergers of gas- 
rich infalling objects at early times, in a number of ways. 
For example, when gas-rich infalling objects merge, the 
gas is likely to form stars with low net angular momen- 
tu m. This process i s refer red to as "dissipative merging" 
by iBekki fc Chlbal (pOOlh and seen (to excess) in early 
cosmologic ally-based galaxy form ation simulations such 
as that of iNavarro et al.l |l995s). Stars formed in this 
way would likely have low metallicity and make up part 
of the inner halo, although the gas might h a ve ex peri- 
enced some pre-enrichment. IBekki fc Chibal (|200lD use 
this process to explain the moderately-flattened, some- 
what metal-enriched inner halo. 

4.4. Dynamical Friction 

There is another way for an accreted satellite to sink 
to the inner halo, even without a radial orbit. Dynam- 
ical friction is a purely gravitational process, caused by 
the formation of a gravitational wake as a massive object 
travels through a density field. It allows satellites, given 
sufficient mass and time before disruption, to transfer en- 
ergy and angular momentum to the dark halo and sink 
toward the galactic center and the disk plane. Dynam- 
ical friction from the halo will permit massive satellites 
which are accreted at a later time to end up in the in- 
ner halo, even if their initial perigalactica are large. The 
efficiency of dynamical friction is a strong function of 
the mass ratio between object and host. In the case of 
a sin gle host, it varies strongly with the satellite mass 



CM 2 ; Binney fc Tremainel [1.987); see also Figure 13 of 



9 Benson (2005) cautions that larger studies are needed to dis- 
entangle the effects of redshift and mass on this quantity. 



I Walker et alT 7l996). so we would expect the more mas 
sive o bjects (which are also likely to form more metal-rich 
stars: Mateo 1998L iTremonti etaHl2004L iLee et al.l[2006L 
lErb et all 120011 to be more successful in reaching the 
solar neighborhood via this process. 

The strong dependence of dynamical friction on the 
mass of the accreted object is important for the inner 
halo: it means that it will be domina ted by debris from 
only a few relatively massive satellites (jHelmi et al.l l2002; 
2003), and m ay well be re l ated to the dumpiness we see 
in E and L. IMeza et al.l (|2005f ) present simulations to 
show that the disruption of a single satellite via dynam- 
ical friction can in fact lead to a number of long-lived 
groupings in energy and angular momentum space. 

Dynamical friction can also affect the spatial distribu- 
tion of the inner halo, as follows. Theory predicts that 
only a few progenitors contributed most of the mass of 
the inner halo, making it less likely that the result of 
these few accretions will be exactly spherical. Accretion 
of hundreds of smaller sub-halos will likely produce a 
spherical outer halo, while it would not be surprising if 
the stars stripped from the few massive objects which 
formed the inner halo produced a somewhat flattened 
inner stellar halo. So the cosmological predictions for 
the growth of the dark halo will quite naturally produce 
a moderately flattened inner halo and a near-spherical 
outer halo. 

4.5. Disk Formation 

As we have noted above, the Galaxy's potential needs 
to settle down (the epoch of major mergers needs to end) 
before we can form either the dynamically cold thin disk 
or our highly flattened halo component, since they would 
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be destroyed by these violent variations in the gravita- 
tional potential. The need for angular momentum con- 
servation in producing the large disks that we see to- 
day also requires this. While the classical picture of 
smooth gas dissipation and infall conserves angular mo- 
mentum, mergers between gas-rich objects at early times 
can transfer both energy and ang ular momentum to the 
dark halo (eg lNavarro et al1l!995f ). Angular momentum 
will not be conserved during this process, which domi- 
nates at early times, unless feedback processes keep gas 
out of the c ores of the merging objects. This was first 
suggested by lNavarro et al.l (|1995l ): we will see below that 
recent simulations have confirmed the result. 

Galaxy formation simulations cannot trace star for- 
mation in detail because the spatial resolution needed 
far exceeds what is computationally possible when 
studying an entire galaxy. Thus average prescrip- 
tions are used to model the complex, multi-phase gas 
physics. In many of the early simulations of galaxy 
formation, the angular momentum of the gas was not 
conserved: gas cooled, collapsed and formed stars 
early, and then mergers tr ansferred its angular mo- 
mentum to t he da rk halo (eg ISteinmetz fc Mullerl Il995l . 
lAbadi et al.l 120031 ) and only small disks were formed 
(the "angular momentum problem" ) . Recent simula- 
tions w hich include a more comple x treatment of feed- 
back (eg lSpringel fc Hernquistll2003h have formed l arger 
more realistic late-type disks ([Robertson et al.l 120061 
IZavala et afll2008l IGovernato et al.ll2007f ). Because the 
accretion of a large satellite will destroy a dynamically 
cold stellar disk, the present-day disk starts to form only 
after major mergers have ceased. However, the star for- 
mation rate in the disks formed in even the most re- 
cent studies tends to peak at early times, soon after 
the last merger (see, for example, lAbadi et al.l I2003L 
IGovernato et alJl2007f ). in contrast to the star formation 
rate in the Milky Way disk, which is th ought to have 
been roughly constant fo r most of its life (|Twarog||1980l . 
iRocha-Pinto et aLll2000t ). More work on modelling feed- 
back and gas flows is clearly needed before the major 
properties of late-type disks are reproduce d. 

In the simulations of lAbadi et al.l (12003) and 
ISpringel fe Hernquist] (|2005l ). the last major mergers be- 
fore disk formation were between gas-rich galaxies. It is 
possible that one or both of the halo's flattened compo- 
nents were created in a way that was intimately related to 
the formation of the disk. In these simulations, the stars 
which were already part of the galaxies in the last ma- 
jor mergers formed a flattened component, and gas from 
these merging objects then settled into a thin disk and 
formed stars. (It is not clear at this stage whether the 
flattened component formed from accreted stars is caused 
by the specific initial conditions of these simulations or 
whether it is a general outcome.) lAbadi et alJ (|2003bl ) 
identify this flattened component (despite its rather large 
scale height of 2.7 kpc) with the thick di sk, which has a 
scale height of ~1 kpc in the Milky Way (Re vle fe Robin! 
120011 iNg et al.l Il997f ) and note that it has kinematics 
"intermediate between halo and disk" : the stars have a 
mean rotation of around 150 km s _1 but a large velocity 
dispersion. This moderately flattened component might 
als o be related to the fla ttened inner halo first suggested 
by iKinman et all (|1965f h As more simulations are pub- 
lished, it will be interesting to see the variation of scale 



height and rotation of the components formed from the 
stars in the satellites which define the orientation of the 
disk. 

It is also possible that the low L±_ group was formed 
at this stage. The gas that came in with them could 
form stars in a flattened configuration that would form 
our low Lj_ group, instead of starting to form the disk. 
Which outcome occurs (rotationally supported disk or 
highly flattened halo) will depend on the angular mo- 
mentum content of this gas. Since, by construction, the 
merger of these satellites was the last major one experi- 
enced by the Galaxy, any flattened stellar configuration 
formed after this time will not be subjected to huge vari- 
ations in the Galaxy's potential, and so it will be able 
to survive. This would explain the fact that the highly 
flattened component has stars which are more metal-rich 
than the overall halo, because stars from the galaxies 
which merged would have already polluted their ISM 
with some metals. While gas on predominantly radial 
orbits in an extremely flattened configuration (like the 
present-day young thin disk) will be likely to "self-cross" 
and shock, forming a more rotationally supported com- 
ponent, the low Lj_ group does not have such extreme 
flattening, so the gas carried in with the satellites which 
experience the last major mergers could still form the 
low L ^ component. 

The accretion and settling to the disk plane of gas and 
subsequent star formation will continue to build the thin 
disk in a process that continues to this day. Orbits of 
inner halo stars will be affected by the potential of this 
highly flattened disk, and could increase the flattening 
of the inner hal o via adiabatic contrac ti on, as originally 
pointe d out by iBinnev fe Mavl (|1986T ). iChiba fe Beers! 
(|2001[ ) have shown that this process is not sufficient to 
form an inner halo with axial ratio as high as 0.7; we 
have seen that the actual flattening of the inner halo is 
closer to 0.6. So adiabatic contraction alone is unlikely 
to form the moderately flattened inner halo. 

4.6. Fashionably late: after disk formation 

After the disk forms, we have another source of dy- 
namical friction which can affect the orbits of objects 
which reach the inner Galaxy. Unlike dynamical friction 
with the halo, this process will be more efficient for ob- 
jects on prograde orbits because the lower relative veloc- 
ity of satellite and disk stars enhances the formation of 
the g ravitational wake and so increases dynamical fric- 
tion ()Walker et all [19961 lAbadi et ail I2003L iMeza et all 
|2005() . We would thus expect to see debris from objects 
which experience dynamical friction with the disk to tend 
to have a net prograde rotation. It is still possible for dy- 
namical friction to bring objects on retrograde orbits into 
the inner Galaxy, b ut the process wil l be s ignificantly 
slower. For example. iBekki fc Freema n (2003) show how 
the angular momenta of stars in a possible progenitor of 
the retrograde, tightly bound globular cluster u> Centauri 
change as it is dragged down into the plane. 

After the merger rate of the early Galaxy settles down 
and the disk forms at around 2=2, the physical condi- 
tions are appropriate for the survival of our flattened 
(low L±) group. As time goes on, the turnaround ra- 
dius continues to grow, and, as noted above, it becomes 
steadily more difficult for objects to reach the inner halo 
without the help of dynamical friction. Satellites on ex- 
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tremely radial orbits can still do so, but will not form 
our low L±_ component because it has quite a range in 
L z . So, dynamical friction is an important part of the 
origin of this component. This fits in well with the result 
shown in Figure [9j more massive objects will also form 
more metal-rich stars on average, so this gives a natural 
explanation of the higher metallicity of this component. 

Can we constrain the number of progenitors whose ac- 
cretion formed this highly flattened component? The 
stars in this component have quite a large range of en- 
ergy and L z , although their Lj_ range is small. The en- 
ergy and angular momenta of stars lost from a satellite 
will be determined by both the orbital energy and an- 
gular momentum of the satellite itself, and the position 
of the stars in the Galaxy when they are stripped. If 
the disruption happens within a few kpc of the Galaxy's 
center, the satellite's finite physical size will ensure that 
stars will gain different amounts of potential energy and 
angular momentum depending on their distance to the 
center when they are stripped from the satellite. Thus 
it is possible that even a single progenitor produced this 
component, if it was disrupted close to the Galaxy's cen- 
ter. It is also possible that more than one satellite was 
disrupted to form it. 

The formation of the disk will flatten the inner 
Galaxy's potential, making it easier for dynamical fric- 
tion with the halo to drag objects down into the plane. 
There are now two options: first that the satellite(s) were 
disrupted before they had a chance to experience dynam- 
ical friction with the disk or thick disk; and second that 
they were disrupted later, in which case dynamical fric- 
tion with both halo and disk would have occurred. If 
only dynamical friction with the halo occurred, then the 
slightly prograde rotation of the highly flattened com- 
ponent is simply a coincidence. This would be easier to 
explain with one or a few progenitors. 

If dynamical friction with the disk itself occurred, we 
have a "built-in" explanation for the small prograde rota- 
tion of this component. Dynamical friction with the disk 
is more efficient for objects on prograde orbits, since they 
will spend more time near the gravitational influence of 
the disk stars. In this case we would expect there to be 
a difference in (L z ) between objects with low L± which 
experienced dynamical friction with the disk, and those 
whose orbits do not spend enough time close to the plane 
to experience disk dynamical friction, and so would be 
found with higher L±. 

At what distance above the plane does dynamical fric- 
tion from the disk become important? We can produce 
a rough estimate for the current epoch by calculating 
where the density of the present-day thin disk and dark 
halo are roughly equal in the solar neighborhood. Using 
a simple isothermal halo with V c — 220 km s™ 1 and r c 
= 2 kpc and assuming a thin disk stellar density of 0.1 
Mq per pc 3 we find that the densities are roughly equal 
at 2.3 thin disk scale heights (z = 800 pc) So we would 
expect stars to need to spend significant time within 1 
kpc of the plane in order to feel dynamical friction from 
the disk. This is quite interesting in the light of the L± 
= 350 result, which is equivalent to a population scale 
height similar to the thick disk's. 

Why should dynamical friction produce a smoother 
distribution for stars close to the plane? The more mas- 
sive satellites would have a larger stellar velocity dis- 



persion, but they will also produce more streams, which 
would need velocity accuracy of order 1 km s _1 to sepa- 
rate. If there was a single progenitor, the range of poten- 
tial energy and angular momentum imparted to stars as 
they were stripped at different galactocentric radii would 
likely produce such a smooth distribution. It is more dif- 
ficult to explain this smoothness in the case of a larger 
number of progenitors. 

Another possibility is that the halo star orbits with low 
L± are close enough to the plane to feel the effect of inho- 
mogeneities in the disk potential due to spiral arms, gi- 
ant molecular clouds and the bar. Scattering and heating 
caused by interactions with these massive objects would 
smooth the velocity distribution of the stars. However, 
the non-circular orbits of these halo stars ensure that the 
stars would not linger long within the gravitational reach 
of these large density inhomogeneities. 

The origin of this low L± group from processes in- 
volving dynamical friction makes the nine globular clus- 
ters with extended horizontal branches and o rbits simi- 
lar to the low L± stars particularly interesting. I Lee et aTl 
(|2007f ) suggest that the clusters with extended horizontal 
branches also have helium-enhanced second-generation 
subpopulations (which are directly observed in several 
cases, notably in u> Centauri). If the clusters were origi- 
nally co res of ancient dwarf galax i es which have since dis - 
rupted (jBekki fc Freeman! I2003L iBekki fc Norrisl l2006ft . 
and the galaxies were sufficiently massive to be affected 
by dynamical friction, then it would not be surprising to 
find them well represented in our low L± group. In fact, 
five out of the six most tightly bound clusters with these 
orbital properties have extended horizontal branches, so 
their orbital properties are consistent with their origins 
as cores of ancient dwarfs. 

In summary, we feel that late accretion, and dynamical 
friction with the dark matter halo (and perhaps also ex- 
isting disk and thick disk stars) is the best way to explain 
the origin of the highly flattened inner halo. 

5. SUMMARY 

Using our sample of metal-poor giants, RHB stars 
and RR Lyraes, we find that the metal-poor halo has 
a highly flattened component, with axial ratio c/a ~ 
0.2, a similar scale height to the thick disk, in the 
neighborh ood of the Sun. This compo nent was sug- 
ges ted by ISommer-Larsen fe Zhenl (Il990f) and predicte d 
bv lPreston et all (|1991l ) and I Green fc Morrison! (|1993f) . 
The highly flattened component has a small prograde 
rotation, but most of its stars are on fairly eccentric 
orbits; this small rotation does not affect its flattening 
significantly. Stars of only moderate metal deficiency 
(—1.5 < [Fe/H] < —1.0) are more likely to be found in 
this component. This component is distinct from both 
the metal weak thick disk, which is primarily rotationally 
supported, and from the moderately flattened inner halo 
detected by many previous workers, which has a much 
less flattened density distribution. 

We suggest that the most likely formation paths for 
this component would be "fashionably late" ones: ei- 
ther the stars formed from the gas accreted in satellites 
whose merger preceded the formation of the disk, or the 
stars were carried into the inner Galaxy by dynamical 
friction acting on their massive progenitor. Other for- 
mation paths such as dissipative mergers at early times 
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could have formed the moderately-flattened halo, but are 
unlikely to have formed our highly flattened component 
because the violent merging activity at these early epochs 
would have destroyed such a coherent component. 

We find that the other stars in our sample, whose or- 
bits reach further from the plane, exhibit quite a clumpy 
distribution in energy and angular momentum, suggest- 
ing that the chaotic conditions under which the early 
halo formed were not violent enough to erase the record 
of the conditions under which they were formed. Any 
truly smooth inner halo is either a very small compo- 
nent or only exists inside the solar neighborhood. This 
clumpy structure also suggests that only a small number 
of objects contributed the majority of stars to the inner 
halo, in line with theoretical expectations. It is possible 
that the progenitor galaxy of u> Centauri was one of these 
objects. 

In addition we find that the RR Lyrae variables in our 
sample are somewhat less tightly bound to the Galaxy, 



on average, than the red giants, and show a more ret- 
rograde rotation. This is consistent with earlier claims 
of a y oung er ou t er hal o such as that of ISearle fc Zinnl 
(119781) andlZinnl (11993) , with the retrograde outer halo 
of iCarollo et al.l (|2007D and with the difference in kine- 
matics bet ween RR Lyrae and B HB stars at the NGP 
detected by iKinman et all (|2007t ). 
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